Direct Evidence of the Discontinuous Character of the Kosterlitz-Thouless Jump 
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It is numerically shown that the discontinuous character of the helicity modulus of the two- 
dimensional XY model at the Kosterlitz-Thouless (KT) transition can be directly related to a 
higher order derivative of the free energy without presuming any a priori knowledge of the nature of 
the transition. It is also suggested that this higher order derivative is of intrinsic interest in that it 
gives an additional characteristics of the KT transition which might be associated with a universal 
number akin to the universal value of the helicity modulus at the critical temperature. 
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I. INTRODUCTION 

The Kosterlitz-Thouless(KT) transition has attracted 
a steady interest since its discovery; ^-^ This is both be- 
cause its unusual characteristic properties and its ap- 
plicability to many systems with two-dimensional (2D) 
character, like, e.g., superfluid/superconducting films. 
The critical properties are given by Kosterlitz's renormal- 
ization equations (KRE)i^ A key characteristic feature is 
the discontinuous universal superfluid jump to zero at 
the transition.^ The generic model for the KT transition 
is the two-dimensional XY model and in this case the 
discontinuous jump is associated with the helicity mod- 
ulusi^i The 2D XY model has been the subject of very 
many computer simulations directed at verifying or dis- 
proving the various characteristics of the KT transition 
as given by KREiSii Although the general consensus is 
that the 2D XY model does indeed undergo a KT tran- 
sition, there have also been claims from simulations that 
the predictions from KRE may not be entirely correct 
In particular, these earlier works have been focused on 
the question whether or not the divergence of the cor- 
relation length (when the transition is approached from 
above), is governed by an essential singularity, as pre- 
dicted by KRE, or by a conventional power-law singular- 
ity. Some evidence for the latter possibility was found 
in Refs. 00- Alternatively, in Ref . m it was argued that 
the divergence is given by an essential singularity which 
is not entirely consistent with KRE. All this reflects the 
difficulty of verifying the precise nature of the transition 
through computer simulations^ 

Similarly, it is difficult to directly verify from simula- 
tions that there is a discontinuous jump in the helicity 
modulus as predicted by KRE. The reason for this is il- 
lustrated in Fig. ^ Simulations can only be performed 
on a finite systems, resulting in the helicity modulus as a 
continuous function without any singularity. As seen in 
Fig- n the drop at the transition gets steeper and steeper 
as the system size is increased. However, the size de- 
pendence is shown to be rather weak. Consequently, al- 
though numerical results like those presented in Fig. 



FIG. 1: Helicity modulus (T) for the 2D XY model for lat- 
tice sizes L = 4 to 64. The transition is signaled by a rapid 
decrease of the helicity modulus in the vicinity of T = 1 as 
temperature is increased. This decrease becomes sharper with 
increasing L. The question is whether or not this rapid de- 
crease develops into a discontinuous jump from a finite pos- 
itive value to zero as predicted by the Kosterlitz-Thouless 
scenario. The data are consistent with such a scenario, but 
on the basis of this type of data alone one cannot rule out 
that the helicity modulus instead goes continuously to zero. 



for the helicity modulus are consistent with a discon- 
tinuous jump for an infinite system size, one cannot on 
this evidence alone, exclude the possibility that the he- 
licity modulus remains continuous in the thermodynamic 
limit. The correctness of the discontinuous character of 
the jump has in practice only been verified from simula- 
tions in more indirect ways which make use of additional 
KRE-predictions for the KT transition like e.g. that the 
leading size dependence of the helicity modulus at the 
critical temperature is logarithmic^ 

Thus one may ask if there exists a more direct way of 
inferring from simulations that there is a discontinuous 
jump, without resorting to more indirect methods which 
use additional KRE predictions. In this paper we show 
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that such a more direct way does indeed exist. 

Our approach is based on the calculation of a higher 
order correlation function. This correlation function ap- 
pears to give an additional and somewhat unexpected 
characteristic feature of the KT transition. We speculate 
that this feature might be associated with a universal 
number akin to the universal value of the helicity modu- 
lus at the critical temperature. 



II. MODEL 

The 2D XY model on a square lattice of the size Lx L 
is defined by the Hamiltonian 



H = cos 

iv) 



L ^ 



(1) 



where J is the coupling strength (set to unity from now 
on), the sum is over nearest neighbor pairs separated by 
the displacement = x or y (we set the lattice spacing 
a = 1), and the phase angle 9i {0 < 9i < 27r) at the lattice 
point i satisfies the periodicity 0i+Lx — Oi+Ly — (^i- For 
generality, we have also included an externally imposed 
global twist across the sample, A = (A^,, Aj^), defined by 
that the summation of the phase difference along the 
X {y) direction equals A^^ (^j/)- The partition function 
Z is given by 



Z 



n 



27r' 



-H/T 



(2) 



and the free energy by F{A.) — —TlnZ. The ground 
state corresponds to the configuration where all spins 
point in the same direction (i.e., all = 0) which means 
that the minimum of the free energy corresponds to 
A = 0. From the symmetry in the Hamiltonian, all odd- 
order derivatives such as 9F/9A|a=o and 9^F/5A'^|a=o 
vanish. Accordingly, in the following we will from stan- 
dard Monte Carlo simulations obtain the two first nonva- 
nishing derivatives of the free energy with respect to A at 
A = 0, i.e., the second order modulus, usually called the 
helicity modulus, (T) = d^F/dA'^\A=o, and the 4th or- 
der modulus (T4} = d'^F/dA'^\/\=o, where (• • •) denotes 
the thermal average. The former (the helicity modulus) 
can be expressed as 

(T) = (e)-:^(.2) 



where 



j2 " ^j)' 



-y 



sm 



and the sum is over all links in one direction. This means 
that (e) is the energy per link and (s) is the average 
current per link in one direction (here taken to be the x 
direction). Note that (s) = by symmetry, in contrast to 
the current-current correlation (s^) which measures the 
current fluctuation. The fourth-order modulus can be 
expressed as 



(L^Ti) = -4(T) +3 



(e)-^((T-(T))2) 



2L^ 



(3) 

In the following, we will measure these two correlation 
functions, (T) and (T4), in standard MC simulations to 
demonstrate the existence of a discontinuous jump of the 
helicity modulus to zero at the phase transition. 



III. STABILITY ARGUMENT 

Since the global minimum of the free energy F{ A) cor- 
responds to zero twist, it follows that -F'(O) < F{A). To 
lowest orders in the A expansion we have 



F(A) = (T) 



A2 



(T4) 



A4 
4! 



This means that (T) > since the lowest order nonvan- 
ishing derivative of the free energy will always dominate 
for small enough A. However, it also implies that the 
next order derivative (T4) likewise has to be > at any 
T where (T) = in the thermodynamic limit. Our ar- 
gument is then the simple observation that (T) cannot 
go continuously to zero at the transition temperature Tc 
if (T4) at the same time approaches a nonzero negative 
value at Tc- But, since (T) is indeed zero in the high- 
temperature phase, this means that if (T4) approaches 
a negative value at Tc then the jump has to be discon- 
tinuous. The point is, as we will show below, that the 
conclusion that (T4) is nonzero and negative at Tc can 
be convincingly drawn from standard MC simulations. 



IV. SIMULATION RESULTS 

Our simulation of the helicity modulus for 2D XY 
model is shown in Fig. ^ up to sizes L — 64. The con- 
clusions which can safely be drawn by analyzing the fi- 
nite size dependence of (T) is that it is finite in the low 
temperature phase and is zero in the high temperature 
phase. Since this is well established we will not discuss 
it further here. As mentioned in Sec. ^ the difficult part 
from a simulation aspect is to determine how the helic- 
ity modulus approaches zero at Tc in the thermodynamic 
limit. 

To this end we will use the stability argument given 
above and instead focus on the size dependence of the cor- 
relation function {Lp'T i) (see Fig. Figure [2Ia) shows 
that (L'^Ti) vanishes in the high-temperature phase and 
goes to a nonzero negative value in the low temperature 
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FIG. 2: (a) The correlation function {L^Ti). As seen this 
correlation function is always negative or zero. In the high-T 
phase it approaches zero from below. In the low-T phase it 
goes to finite negative values, as is apparent from the data 
points somewhat below T = 1. The most interesting feature 
is the divergent dip in the vicinity of T = 1. This divergence 
at Tc is the manifestation of the phase transition for this cor- 
relation function. The inset illustrates that the half-width of 
the dip goes to zero for large L. (b) The 4th order modulus 
(T4) for lattices sizes L = 12 to 64. As seen this modulus 
goes to zero both below and above Tc. However, the crucial 
point is that the depth of the dip remains finite, as is shown in 
the inset where the data are plotted against 1/L. The linear 
extrapolation to L = 00 gives (T4) ~ —0.130 ± 0.005 at Tc 



phase as the system size becomes larger. The interesting 
features is the divergent dip. This occurs in the region 
where the helicity modulus for finite L goes rapidly to- 
wards zero. Thus this singularity in (L'^T^) can safely be 
associated with Tc. The conclusion from Fig.jS^a) is that 



(L^T4) diverges at Tc for L ^ 00, goes to zero above, and 
goes to a negative nonzero value below. This conclusion 
is further supported by the fact that the half-width of the 
divergent dip decreases towards zero for L — > 00 as illus- 
trated in the inset of Fig. [5fa). The crucial point is now 
what this divergence in the correlation function (L'^T-j) 
implies for (T4). This is shown in Fig.|2Ib), which gives 
the result for L > 12. The inset in Fig. |2tb) gives the 
depth of the dip as a function of size. As is apparent from 
the inset, the depth goes to a finite value in the thermo- 
dynamic limit. Linear extrapolation in 1/L to L = 00 
gives the value 0.130 ± 0.005. Thus we conclude that 
(T4) is indeed negative and nonzero precisely at Tc. Us- 
ing the stability argument in Sec. lIIII this means that the 
helicity modulus has to be positive and nonzero precisely 
at Tc. However, above Tc the helicity modulus is zero 
from which follows that the helicity modulus has to be 
discontinuous at Tc. 

One may also note from Fig. I3b) that the position of 
the minimum of the dip decreases towards lower T with 
increasing L. However, the approach towards the known 
value of Tc (~ 0.89) is rather slow making a precise de- 
termination of Tc based on (T4)-data less advantageous. 



V. FINAL REMARKS 

We have shown that the discontinuous character of the 
jump of the helicity modulus at the transition for the 2D 
XY model can be established from MC simulations of the 
4th order modulus given by Eq.lQ. It seems very likely 
that the singularity found here for the 4th order modulus 
in case of the 2D XY model is part of the general KT 
transition scenario. Thus we suggest that it sometimes 
may be advantageous to study the singularity of this 4th 
order modulus in simulations, when trying to determine 
if a transition for a particular model is of KT type. It is 
also tempting to speculate that the finite value of (T4) w 
—0.130 at Tc is associated with a universal number, akin 
to the universal value of of the helicity modulus (T) = 
2Tc/^. 
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